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Abstract 

We evaluate both the massless and the massive Dirac quasi-normal mode frequencies in the 
Schwarzschild black hole spacetime using the WKB approximation. For the massless case, we find 
that, similar to those for the integral spin fields, the real parts of the frequencies increase with the 
angular momentum number k, while the imaginary parts or the dampings increase with the mode 
number n for fixed n. For the massive case, the oscillation frequencies increase with the mass m of 
the field, while the dampings decrease. Fields with higher masses will therefore decay more slowly. 
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I. INTRODUCTION 

In a series of papers, Finster and collaborators studied the 

interaction of the Dirac field with gravity. While they have found stable particlelike solu- 
tions in the Einstein-Dirac-Maxwell system hi yf , they have also proved the non-existence 

nn 

of time-periodic solutions in various black hole spacetimes jj, |5||. This is quite surprising 
because it means that Dirac particles, including electrons and neutrinos, cannot remain on 
a periodic orbit around a black hole. In J3|, they also showed that there are no spheri- 
cally symmetric black hole solutions in the Einstein-Dirac-Maxwell system other than the 
Reissner-Nordstrom one. This indicates that if a "cloud" of Dirac particles collapses grav- 
itationally, the particles must eventually vanish inside the event horizon of a black hole 
or escape to infinity. It is therefore interesting to see how Dirac fields evolve in curved 
background spacetimes. To understand this, Finster and collaborators considered the mas- 
sive Dirac field in the general Newman-Kerr black hole spacetime and derived the so-called 



power-law tail behavior of the field at very late time in 

The evolution of wave fields around a black hole consists roughly of three stages fu . 
The first one is an initial wave burst coming directly from the source and is dependent on 
the initial form of the original wave field. The second one involves the damped oscillations 
called the quasi-normal modes, which do not depend on the initial values of the wave but are 
characteristic of the background black hole spacetimes. The normal mode frequencies have 
complex values because of radiation damping. The last stage is the power-law tail behavior 
mentioned above. It arises because of the backscattering by the long-range gravitational 
field. 

Here we would like to concentrate on the intermediate stage of the evolution of the 
Dirac field where the quasi-normal modes dominate. As a first step we would calculate 
the massive Dirac quasi-normal mode frequencies in the Schwarzschild black hole spacetime 

n 

in this paper. The study of black hole quasi-normal modes has a long history 12]. It 
started with the black hole stability problem which concerned the evolutions of black hole 
perturbations. Since these perturbations are effectively wave fields of spins 0, 1, and 2, the 
evaluations of the quasi-normal mode frequencies are mainly done for these cases. As far 
as the Dirac field is concerned, we are not aware of any work done in this direction. So we 
hope that our calculation can serve to fill in this gap. 
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There are quite a few methods in evaluating the quasi-normal mode frequencies. Most of 
them are numerical in nature. However, a powerful analytical WKB scheme was devised by 
Schutz and Will and was subsequently extended to higher orders in 14 1. Comparing 



with other apples this S cW , show. t o he _ for the lowing .odes 
Since we are interested in the evolutions of the Dirac fields in the intermediate stage, which 
is dominated by the low-lying modes, this WKB scheme is good enough for our purpose. 

In the next section, we consider the Dirac equation in the Schwarzschild spacetime and 
its reduction into a set of Schrodinger-like equations with a particular effective potential. 
In Section III, we evaluate the quasi-normal frequencies using the WKB scheme for the 
massless case. The massive case is considered in Section IV. Conclusions and discussions 
are presented in Section V. 

II. DIRAC EQUATION IN THE SCHWARZSCHILD SPACETIME 

Consider the Dirac equation in a general background spacetime 

[ 7 a e/(<9 M + rg+m]^ = 0, (1) 

where m is the mass of the Dirac field, and e a ^ is the inverse of the tetrad e^ a defined by 
the metric g^ v , 

9[tv = Vab&ii &v i (2) 

with rjab = diag(— 1, 1, 1, 1) being the Minkowski metric. 7 a are the Dirac matrices 




(3) 

where a 1 are the Pauli matrices. is the spin connection given by 

T, = \b\l b We bv ,^ (4) 

where e&„ ;/Lt = d^thv — ^ u e ba is the covariant derivative of with r" being the Christoffel 
symbols. 

In the Schwarzschild spacetime, 

^ = _A^^ + / 1 _^y 1 dr 2 + r 2 M 2 + r 2 sin 2 9 #2) (5) 
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, r, r sin# . 



where M is the mass of the black hole. Here we can take the tetrad to be 

// 2M\ 1/2 / 2M\" 1/2 \ 
e/ = diag^l-— ) ,(!-—) 

Then the spin connection defined by Eq. (jlj) is 

i / 2M\ 1/2 /l M 

7Vr, = V i— ' 



(6) 



r 2r(r-2M) 



2r 



cote. 



(7) 



Therefore, the Dirac equation becomes 



/ _ 2M_y 1/2 d_ + 1 /_ 2M\ 1/2 (9_ I 
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+7' 



rsin# 



.1/ 



2r(r - 2M) 
* = 0. 



The equation can be simplified by defining 



2M\- 1/4 
r / 



(8) 



(9) 



Then Eq. (JHJ) becomes 



V r / at \ r J \or r J 



' d l 

r) V 00 + 2 



cot0 + 7' 



/ 1 



\rsin# 







+ m 



$ = 0, 



(10) 



which is closely related to the Dirac equation in flat spacetime with a central potential Q| 
We can try the ansatz, 

' iG<±)(r) (±)//j ,s 



-iEt 



where for j — I + 1/2, 



and for j = I — 1/2, 



/ / j+!72+m v m-l/2 
V 2Z+1 r ' 
/ l+l/2^m ym+l/2 

V V 2/+1 r * 



(11) 



2Z+1 r ' 



l+l/2+mym+l/2 



(12) 



(13) 



Then the radial equations for and can be simplified to 



d 



i?(±) 



2M 



i?(±) 



m —Kr±-)/r ) \ G^> 
4 



-£ 
£ 



/r(±) 



(14) 



where we have made a coordinate change 



r + 2Mln 



2M 



- 1 



«(±) 



and the constant 

' -(j + 1/2), j = Z + l/2 

a +i/2), j =1-1/2 

Note that and are negative and positive integers, respectively. 

To further simplify the radial equations, we consider separately the (+) and 
First, for (+), we make a change of variables ^| 

/ 



( FW 



where 



sin(0 (+) /2) cos(0 {+) /2) 
^cos(0 (+) /2) -sin(0 (+) /2) 

(+) = tan" 1 (mr/ 1 « (+) |). 



Then Eq. (jTljl becomes 



d 

dr* 



#(+) 



1 



2M //« 



+ m 2 



\ 2E\ r J k 2 ( ) + m 2 r 2 I I x 



i?(+) 

b (+) "r "' ' / \ ! u / \ G (r 
This equation can be simplified further if we make yet another change of variable, 



1 
-1 

-1 




G« 



-i 



and then 



where 



_d_ 

dr* 



r * + 2E tan 



G« G(+) 



mr 



c (+)l 



G (+) 
__p(+) 



1 ZiV'i 



(+) 



2MJ im 



+ m 2 



1 + 



"(+) 



These equations for and G^ can be decoupled giving 



dfl 



where 



V W = ± d -^ + W? +y (25) 



dr. (+) 

For the case of (— ), similar procedures can be carried out giving 



'4 + %)^ = ^ ( 26 ) 



df 2 



+ V H2 )G^ = E 2 &~\ (27) 



where 



with 



V W = ±^ + Wl_ ) , (28) 



m 2 



W H = — ! (29) 

1 L fl _ 2AA [ „ mK (~) 1 

1 2£ V r y ^^+raV 2 y 

From the forms of W(±) in Eqs. (|22[1 and ([29)1 . we see that the (+) and (— ) cases can be put 
together giving 

(-i + v >) p = E2p ' (30) 

Vi : 2 = ±^ + M /2 , (32) 



where 



with 



_ ^JMzj +m ' 2 

i , 1_ / 1 _ 2M \ / 

1 " l " 2£ V r J U 2 +m 2 r : 

Al/2( /t 2 +m 2 r 2 ) 3/2 



(33) 



r 2 (/t 2 + m?r 2 ) + m«;A/2i?' 
where A = r(r — 2M). Here k goes over all positive and negative integers. Positive integers 
represent the (+) case with 

k = j + - and j = l + - (34) 
while negative integers represent the (— ) case with 

/-,' = - I j + ^ and j = I - ~ (35) 
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From the Schrodinger-like equations in Eqs. (|3U|) and pip, we shall evaluate the corre- 
sponding quasi-normal mode frequencies. Note that V\ and V2, which are related as shown 
in Eq. (|32p. are supersymmetric partners derived from the same superpotential W Q|. It 
has been shown th at p otentials related in this way possess the same spectra of quasi-normal 
mode frequencies |20|. Physically this just indicates that Dirac particles and antiparticles 
have the same quasi-normal mode spectra in the Schwarzshild black hole spacetime which 
is quite reasonable here. We shall therefore concentrate just on Eq. (j3*Uj) with potential V\ 
in evaluating the quasi-normal mode frequencies in the next sections. 



III. QUASI-NORMAL MODE FREQUENCIES FOR THE MASSLESS DIRAC 
FIELD 

In this section we shall evaluate the quasi-normal frequencies for the massless Dirac field 
using the WKB approximation. For the massless case, the radial equation (Eq. (jHOj)) is 
simplified to 

/ d 2 \ ~ 

+ V(r, k) )F = E 2 F, (36) 



\ dr 2 1 
where 

V(r,«) = J^^[|«|A 1 / a -(r-3)] > (37) 

r* = r + 2 ln(r/2 — 1), and A = r(r — 2). Note that we have written V\ as V because we shall 
not work with V2, which will give the same spectrum of quasi-normal mode frequencies, as 
discussed in the last section. Also we have used the mass M of the black hole as a unit of 
mass and length to simplify the notation. 

The effective potential V(r, k), which depends only on the absolute value of k, is in the 
form of a barrier. As shown in Fig. ^ the peak of the barrier gets higher and higher as |ac| 
increases. From Fig. EJ we see that the position of the peak 

r max {\K\ -> 00) -> 3, (38) 

and the maximum values of the potential barrier increase indefinitely with \k\. In fact, 

K 2 

V(\k\ -> oo,r max ) = — (39) 



as can be easily derived from Eq. ()37p. 
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FIG. 2: Variation of the positions of the peaks (left) and the peak values (right) of the effective 
potential for the massless Dirac field with 

To evaluate the quasi-normal mode frequencies, we adopt the WKB approximation devel- 
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TABLE I: Massless Dirac quasi-normal mode frequencies. 



\k\ 


n 


Re(E) 


lm(E) 


1 








176 


-0 


100 


2 








379 


-0 


097 




1 





354 


-0 


299 


3 








574 


-0 


096 




1 





556 


-0 


293 




2 





527 


-0 


497 


4 








767 


-0 


096 




1 





754 


-0 


291 




2 


o 


'von 

730 


-U 


491 




3 


o 


700 


-0 


696 


5 








960 


-0 


096 




1 





950 


-0 


290 




2 





930 


-0 


488 




3 





904 


-0 


689 




4 





872 


-0 


894 



oped by Schutz, Will, and Iyer 



n 



14|,|lj|. This analytic method has been used extensively 
in various black hole cases. Comparing with other numerical results, this method has been 
found to be accurate up to around 1% for both the real and the imaginary parts of the 
frequencies for low-lying modes with n < I, where n is the mode number and I is the angular 
momentum quantum number. The formula for the complex quasi-normal mode frequencies 
E in the WKB approximation, carried to third order beyond the eikonal approximation, is 
given by ^| 

(40) 



E 2 = [V + (-2<) 1/2 A] - i(n + i)(-2V ") 1 / 2 (l + fi), 



where 



A 



n 



i 



1 Vo 



(4)' 



-2V0') 1 / 2 }8\Vo)U + a ) 288 V V 



1 



Vr 



-2V ") 1 6912 \V C 



(77 + 188a 2 




(41) 



(51 + 100a z 
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FIG. 3: Massless Dirac quasi-normal mode frequencies 



v {4) V i fv'"v {5y 



( 67 + 68 a 2 ) + ~~ T77^ (19 + 28a' 



2304 \ VS ' v 288 I ^ 



Here 



a = n + — , n=< (43) 




= — ■ (ID 

Plugging the effective potential in Eq. (J37|) into the formula above, we obtain the complex 
quasi-normal mode frequencies for the massless Dirac field. The values for < n < k are 
listed in Table The values for negative n are related to those with positive n by a reflection 
of the imaginary axis so we do not list them out here. The quasi-normal mode frequencies 
for positive k are plotted in Fig. EJ The variations of the mode frequencies here are similar to 
those with integral spin fields in the Schwarzschild black hole spacetime The real part 
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Ke(E) of the frequencies decreases as the mode number n increases for the same angular 
momentum number k. On the other hand, the magitude of the imaginary part lm(E) of the 
frequencies increases with n. This indicates that the quasi-normal modes with higher mode 
numbers decay faster than the low-lying ones. Hence, the low-lying quasi-normal modes, 
possibly with different values of \k\, are most relevant to the description of the evolution of 
a massless Dirac field in the Schwarzschild black hole spacetime. 



IV. QUASI-NORMAL MODE FREQUENCIES FOR THE MASSIVE DIRAC 
FIELD 



In this section we shall evaluate the quasi-normal mode frequencies of the massive Dirac 
field. Again we start with the radial equation (Eq. ([30])). 

/ d 2 



df- 



+ V(r, K,m,E)\F = E F, 



(45) 



where f* = r + 2 m(r/2 — 1) + (l/2£')tan 1 (mr/K), and 



V(r, k, m, E) 

AVy + mV) 3 / 2 
(r 2 (ft 2 + m 2 r 2 ) + mKA/2E) 2 

A3/2( K 2 + m 2 r2) 5/2 



A 1/2 (k 2 + m 2 r 2 f /2 + ((r - l)(/t 2 + m 2 r 2 ) + 3m 2 rA) 
2t(k 2 + m 2 r 2 ) + 2m 2 r 3 + rriK(r - 1) /E 



_ ,. , , , ,, , . . (46) 

(r 2 (K 2 + m 2 r 2 ) + mKA/2E) 3 l J 

Note that the effective potential depends not only on m but also on the energy E. This will 
complicate our evaluation of the quasi-normal mode frequencies using the WKB approxima- 
tion formula in Eq. ()40|) because there are E dependence on both sides of the equation. 



A. Properties of the effective potential 

Here we analyse the dependence of the effective potential on the parameters m, k, and 
E. First, its dependence on m is showed in Fig. 0] with energy E = 1 and with k = 1, 2, — 1 
and —2. For small values of m, the potential is still in the form of a barrier, but with the 
asymptotic value 

V(r -> oo) = m 2 . (47) 

As m is increased, the peak of the potential also increases but does so very slowly. Eventually, 
the height of the peak is lower than the asymptotic value m 2 . In this case there will be no 
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FIG. 4: Variation of the effective potential for the Dirac field with mass m for k = 1,2, —1 and 
-2. 

quasi-normal modes because for tunneling to occur E 2 of the Dirac field must be smaller 
than the peak value of the potential. This cannot happen since the energy of the Dirac 
field is always larger than its mass m. When m is increased further, the peak disappears all 
together, and the potential barrier turns effectively into a potential step. Note that there is 
a caveat here. The energy E is unknown at the start of the calculation (here we have taken 
a typical value E = 1) and has to be determined self-consistently. Hence, Fig. 0] (similarly 
for the other figures below) should only be taken as indicative of the general behaviors of 
the effective potential V. 

The effective potential for the massive Dirac field also depends on the energy E, and we 
show its dependence in Fig. ED E is varied from its minimum value m. For positive k, the 
value of the potential increases with E, while for negative k, it decreases. The variations 
become smaller and smaller, and we can see from the figure that the potentials are nearly 
indistingusable for E = 5 and E = 100 or higher. The general behaviors of the potential 
remain the same for all these values of E. Indeed, from the form of the potential in Eq. (|46j). 
E appears all in the denominators. The terms involving E can never get large enough to 
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FIG. 5: Variation of the effective potential for the Dirac field with energy E for m = 0.2 and 0.4, 
and k = 1 and —1. 

change the general behaviors of the potential since E cannot be smaller than m. 

From the dependences of the potential with m and E, we can also estimate the maximum 
values of m above which quasi-normal modes cannot occur. From the discussions above, we 
see that quasi-normal modes exist only when (i) the peak value of the potential V(r = r max ) 
is larger than m 2 , and (ii) E 2 of the field is smaller than this peak value. Therefore, one can 
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TABLE II: Estimation of the maximum values of the mass m and [i (= m/ k) of the Dirac field 
above which quasi-normal modes cannot occur. 



«| m max (positive k) [i m ax (positive k) m max (negative k) Umax (negative k) 

1 0.224 0.224 0.333 -0.333 

2 0.453 0.226 0.572 -0.286 

3 0.696 0.232 0.819 -0.273 

4 0.944 0.236 1.067 -0.267 

5 1.192 0.238 1.316 -0.263 




FIG. 6: Variation of the effective potential for the Dirac field with k for m = 0.2 and 0.6. 
estimate this maximum value m max from the relation 

which can be solved numerically. The result is tabulated in Table |Hj where we have also 
given the maximum values of \x = m/n which will be useful in the calculation below. 

Finally, in Fig. |H1 we show the dependence of the effective potential with the angular 
momentum quantum number k. We see that as k increases, the behaviors of the potential 
approach to that of the massless one, as shown in Fig. regardless the mass of the field. 
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This can be understood by writing the potential in Eq. (|4l)j) as 

V(r, K, m = Kfj,, E) 

-w^MrSmf I 2r(1 + " V) + 2 " v + Kr " • (49) 

where \i is as defined above. The dependence of mass is all through \i. As |/c| — > oo, \i — ► 
regardless the value of m, and the effective potential will approach to that of the massless 
one. Also from Table ITT1 the maximum magnitudes of fi are around 0.2 ~ 0.3 regardless the 
value of k. n can thus be used as a small parameter for expansions as we shall do below. 



B. Evaluation of the quasi-normal mode frequencies 

The evaluation of the massive Dirac quasi-normal mode frequencies is more complicated 
than the massless ones due to the presence of the energy dependence on both sides of the 
Eq. ()40j) . Since the maximum magnitudes of /i are only around 0.2 ~ 0.3, we can try to 
obtain the quasi-normal mode frequencies in the WKB approximation as power series of \i 



2l| and 



221. 



for given values of k similar to the power series expansions in 

We first express the position of the peak of the effective potential as a series up to order 

rmax = r + mi + r 2 /i 2 + r 3 /i 3 + r 4 /i 4 + r 5 /i 5 + r 6 /i 6 

= r + E, (50) 



by requiring 



= V frv 



= V'(r ) + Ey"(r ) + isV"(r ) + ^ 3 V^(r ) 

2 o 

+ ^£V 5 )(r ) + Y^V^ro) + 4 SV(7)(ro) " (51) 
Since vq is just the position of the peak for the massless case obtained in Section UTT] we can 
evaluate the coefficients rVs order by order from this equation. For example, for \k\ = 1, we 
obtain 

/0.182N , / 0.0755\ „ /1.07 0.0322X 

W = 2.42 - /, (_) + ^ (3.66 + — ) - ^ (_ + _) 
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+// 4 (18.8 - 
+/i 6 (ll7 



0.252 0.0142\ 



E 2 
8.16 



E 4 ') r " V E 
0.597 0.00293^ 



A 4 



0.526 



0.624 0.00638 
+ 



£ 3 



£ 5 



(52) 



£ 2 £ 4 £ 6 

The coefficients in the above expansion involve the unknown E. We can likewise expand E, 



E = E + E 1 fi + E 2 fi 2 + E 3 fi 3 + £ 4 /i 4 + £ 5 /i + £ 6 /i 6 . 



(53) 



We plug this expansion back into Eq. ()52jl. and expand again up to order /i 6 . Using this 
new expansion for r maX} we can then expand the derivative of the potential Vq . Note that 
the derivatives are performed with respect to f*, and 

1 n 



_d_ 

dr* 



1 - 2/r + 2E(1 + /i 2 r 2 ) 



d_ 

dr 



2E 



AE 2 



1 - 



~(l-;) (4-4r + r>-4EVV 

_ ^ 3 ( 4 _ 4r + r 2_ 8jB 2 r 4 )At , 

rj 

. (16 - 32r + 24r 2 - 8r 3 



16£ 4 r 2 
1 



+ (1 - A8E' 2 )r 4 + A8E 2 r 5 - 12£ 2 r 6 + 16-EV)// 



2^6 



1 /. 2^ 3 



64£ 6 r 4 



1 - 



(16 - 32r + 24r 2 - 8r 3 



+ (1 - 64£ 2 )r 4 + 64E7V - lQE 2 r b + 48EV)/i b + 



dr 



(54) 



We put all these expansions back to Eq. (J40j) and solve the coefficients E^s self-consistently 
order by order in \x. The results are tabulated in Tables ITTT1 and ITVl 

With these coefficients we can evaluate the quasi-normal mode frequencies by putting in 
the value of \x for both positive and negative values of k. From Table ITT1 the maximum value 
of \i that could have quasi-normal modes is 0.224 for k = 1 and increases to 0.238 for k — 5. 
In fact, the value will increase for larger values of k, but converges to the value ~ 0.25. 
Similarly for the variation of fi m ax with negative k, the magnitude decreases from 0.333 for 
k = —1 to 0.263 for k = —5, and will converge to ~ 0.25 for large values of k. Therefore, 
we shall try to evaluate the frequencies for from to 0.2. 

First we examine the frequencies for k = 1 and n = 0. For m = 0, we have 



E = E n 



0.176-0.1001 



(55) 
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TABLE III: Real parts of the coefficients of the expansions in powers of fj, for the quasi-normal 
mode frequencies with |re| = 1 to 5. 



n\ 


n 


Re(E ) 


Re(£?i) 


Re(E 2 ) 


Re(£ 3 ) 


Re{E 4 ) 


Re(£ 5 ) 


Re(£ 6 ) 


1 





0.176 


-0. 


168 


0.521 


-0.393 


-0.00367 


-3.13 


74.2 


2 





0.379 


-0. 


166 


1.49 


0.158 


1.01 


-0.728 


-0.845 




1 


0.354 


-0. 


136 


0.478 


-0.510 


-3.81 


9.23 


-24.0 


3 





0.574 


-0. 


166 


2.43 


0.335 


1.84 


0.721 


2.75 




1 


0.556 


-0. 


152 


1.55 


-0.366 


-1.92 


-3.30 


-8.91 




2 


0.527 


-0. 


129 


0.463 


-0.455 


-5.84 


16.6 


-93.0 


4 





0.767 


-0. 


166 


3.34 


0.404 


2.56 


1.13 


6.85 




1 


0.754 


-0. 


158 


2.61 


-0.115 


0.291 


-2.31 


-2.43 




2 


0.730 


-0. 


144 


1.50 


-0.494 


-5.64 


-0.255 


-49.8 




3 


0.700 


-0. 


126 


0.460 


-0.427 


-7.39 


20.8 


-152 


5 





0.960 


-0. 


166 


4.23 


0.437 


3.24 


1.28 


9.98 




1 


0.950 


-0. 


161 


3.62 


0.0610 


1.80 


-0.952 


2.32 




2 


0.930 


-0. 


152 


2.60 


-0.358 


-3.02 


-3.28 


-23.4 




3 


0.904 


-0. 


138 


1.46 


-0.514 


-8.82 


4.03 


-103 




4 


0.872 


-0. 


124 


0.462 


-0.408 


-8.81 


23.5 


-204 



which is the massless case we have obtained in the last section. 
For m — 0.1, we have 



Eq — 


0.176 -O.lOOi, 


E lf i = 


-0.0168 -0.007112, 


E 2f i 2 = 


0.00521 + 0.00452i, 


E 3f i 3 = 


-0.000393 + 0.000375i, 


E,fi 4 = 


-0.00000367 + 0.000138i 


E 5f i 5 = 


-0.0000313 - 0.0000859i, 


£ 6/ i 6 = 


0.0000742 - 0.0000140i, 
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TABLE IV: Imaginary parts of the coefficients of the expansions in powers of \x for the quasi-normal 
mode frequencies with |«| = 1 to 5. 



k\ 


n 


lm{E ) 


Im(Ei) 


lm(E 2 ) 


lm{E 3 ) 


lm(E 4 ) 


Im(£ 5 ) 


lm(E e ) 


1 





-0.100 


-0.0711 


0.452 


0.375 


1.38 


-8.59 


-14.0 


2 





-0.0965 


-0.0282 


0.653 


0.480 


1.38 


1.59 


-3.47 




1 


-0.299 


-0.0771 


1.16 


-0.0876 


-5.43 


-16.9 


7.54 


3 





-0.0963 


-0.0185 


0.690 


0.368 


1.21 


1.46 


4.76 




1 


-0.293 


-0.0533 


1.64 


0.391 


1.38 


-4.30 


22.6 




2 


-0.497 


-0.0770 


1.72 


-0.172 


-13.3 


-18.9 


37.7 


4 





-0.0963 


-0.0139 


0.704 


0.291 


1.13 


1.14 


6.02 




1 


-0.291 


-0.0406 


1.85 


0.516 


2.87 


0.442 


21.0 




2 


-0.491 


-0.0627 


2.38 


0.186 


-2.21 


-11.0 


59.4 




3 


-0.696 


-0.0771 


2.27 


-0.204 


-20.5 


-19.6 


52.6 


5 





-0.0963 


-0.0111 


0.710 


0.239 


1.08 


0.914 


6.23 




1 


-0.290 


-0.0327 


1.96 


0.517 


3.20 


1.60 


19.9 




2 


-0.488 


-0.0521 


2.76 


0.399 


2.25 


-4.00 


53.6 




3 


-0.689 


-0.0672 


3.01 


0.0573 


-7.44 


-15.0 


81.5 




4 


-0.894 


-0.0773 


2.80 


-0.223 


-27.5 


-19.8 


53.6 



and 

E = 0.164 — 0.102* (57) 

Hence, if we take three significant figures, both the real and the imaginary parts of Es^i 3 
and higher order terms do not contribute to the frequency. The result in Eq. (JH7j) is thus 
accurate up to at least three significant figures. 
For m = 0.2, we have 

E = 0.176 -O.lOOi, 
E x n = -0.0337- 0.01422, 
E 2 fx 2 = 0.0208 + 0.01811, 
E 3l x 3 = -0.00315 + 0.00300i, 
E^ 4 = -0.00000587 + 0.00221i 
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FIG. 7: Massive Dirac quasi-normal mode frequencies 
£ 5/ i 5 = -0.00100 - 0.00275i, 

E 6 fi 6 = 0.00475 - 0.000893i, (58) 

and 

E = 0.164 — 0.0947i (59) 

Here we see that even the highest order term E 6 fx 6 will contribute to both the real and 
the imaginary parts of E. If we require our result to be accurate up to three significant 
figures, we need to go to orders higher than /A However, we have only expanded terms up 
to we are thus forced to disregard this result for m = 0.2, even though there should be 
a quasi-normal mode around that value. 



19 



In summary, for different values of k and n, we increase the values of m (or /i) from 
up to the point that three significant figure accuracy cannot be maintained with terms up 
to only /i 6 , and then we stop at that value. The results for k = 1 to 5 are plotted in Fig. [7| 
Comparing with Table HT1 we see that the largest value of m for each k and n that we present 
in Fig. [7|is in fact very close to the maximum value m max . On the other hand, the trends 
of the variations of E with m are similar for negative k so we do not plot them out. 

As we can see from Fig. [?[ the real parts of the quasi-normal mode frequencies in general 
increase with the mass of the field. This can be understood from the fact that when the 
mass of the field increases, the peak of the effective potential gets higher. The real part 
of the energy should therefore increase. We also see that the magnitudes of the imaginary 
parts decrease with the mass. We thus expect Dirac fields with larger masses (same k and 
n) to decay slower during their evolutions. 



V. CONCLUSIONS AND DISCUSSIONS 



We have evaluated both the massless and the massive Dirac quasi-normal mode frequen- 
cies in the Schwarzschild black hole spacetime using the WKB approximation. Only the cases 
with k > n, where k is the total angular momentum number and n is the mode number, are 
considered since these are the cases known to be accurate for the WKB approximation with 

n 

the integral spin fields [15]. For the massless Dirac field, we see that the real parts of the 
quasi-normal mode frequencies increase with k, but they decrease with n for fixed k. As for 
the imaginary parts, the magnitudes of which increase with n. Hence, the fields with n = 
will decay most slowly. 

For the massive cases, the evaluation of the quasi-normal mode frequencies is complicated 
by the fact that the effective potential also depends on the energy E. A further approxima- 
tion is needed to carry through the calculation in which we make perturbative expansions 
for all the relevant quantities in powers of the parameter fi = m/n where m is the mass of 
the field. In this way we are able to obtain the numerical values of both the real and the 
imaginary parts of the mode frequencies up to three significant figures. In general, when the 
mass of the field is increased, the real parts of the frequencies increase while the magnitudes 
of the imaginary parts decrease. Fields with higher masses will therefore decay slower. This 
result is similar to that of the massive scalar field 
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The calculation carried out here can be extended to other black hole spacetimes, for ex- 
ample, the charged Reissner-Nordstrom black hole or the rotating Kerr black hole. However, 
for the Kerr black hole, the situation is further complicated by the fact that there is no close 
form for the effective potential. Furthermore, we can consider dilaton black holes or black 
holes in the anti-de Sitter space both of which are closely related to string theories. 

As we have mentioned ealier, our main concern is about the evolution of Dirac fields 
in curved spacetimes, especially black hole spacetimes. After calculating the quasi-normal 
mode frequencies, our next step would be to understand the late-time tail behaviors of the 
fields. Together with the Green's function method, we should be able to get a more complete 
picture of the evolution, and hopefully towards a better understanding of the gravitational 
collapse of Dirac (or fermionic) matter. 
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